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ABSTRACT 

We describe a generalization of the asymptotic calculation of the tidal torques ex¬ 
perienced by a massive star as a result of a companion in circular orbit originally 
considered by Zahn (1975,1977) to the case of a rotating star when the forcing fre¬ 
quency is small and in the inertial regime, that is it is less than twice the rotation 
frequency in magnitude. The results confirm the presence of a strong toroidal mode 
resonance feature for retrograde forcing and also, with a simple description of the con¬ 
vective core, the presence of some core inertial mode features in the response. These 
were found numerically by Savonije and Papaloizou ( astro-ph/9706186 ). 
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1 INTRODUCTION 

In two recent papers, ( Savonije, Papaloizou and Alberts, 1995: SPA, Savonije and Papaloizou, 1997: SP) we studied the 
response of a uniformly rotating massive star to the tidal forcing due to a companion in a circular orbit (e.g. a massive X- 
ray binary). To make the problem tractable, centrifugal forces were neglected, allowing a spherically symmetric equilibrium, 
but Coriolis forces were retained allowing normal modes governed by rotation to enter into the response. This problem is 
complex because the dense spectrum of normal modes causes excitation of short wavelength responses which lead to numerical 
difficulties. As a result the analysis of SPA could not be extended into the inertial regime in which the forcing frequency is 
less than twice the rotation frequency in magnitude. In SP numerical results appropriate to the inertial regime for modest 
stellar rotation rates have been obtained through the introduction of a viscosity which provides numerical damping of the 
shortest wavelengths in the response in convective regions. 

The first calculations of the tidal torques experienced by a non-rotating massive star as a result of a companion in a circular 
orbit were performed by Zahn (1975,1977) who developed an asymptotic approach valid for low forcing frequencies. Papaloizou 
and Savonije (1984, 1985) found that this could give reasonable results provided the star was unevolved with the convective 
core not being too small in size. 

In this paper we investigate a generalization of the asymptotic calculation of the tidal torques to the case of a rotating 
star when the forcing frequency is small in magnitude and in the inertial regime. Such an analysis is of interest in order to 
complement the numerical work which runs into difficulties because of the short wavelength nature of the response at low 
frequencies. We find that some of the phenomena found by SP, in particular the strong toroidal mode resonance, is present in 
the asymptotic analysis as well and persists in the low frequency limit. This may be important for retrograde forcing leading 
to a more rapid synchronization of spin and orbit than would be expected from consideration of non-rotating stars. 

In section 2 we give the basic equations for a uniformly rotating star subject to tidal perturbation due to a companion in 
circular orbit, formulating a single equation giving the response of the adiabatic interior. In section 3 we consider the response 
of the convective core, indicating how inertial modes may be excited there. We also consider the radiative exterior indicating 
how use of the traditional approximation leads to a separable problem for excited generalized g and r modes in that region. In 
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section 4 we consider the WKB approximation for the radiative region with outgoing wave conditions and how toroidal mode 
resonances lead to a global response. We go on to evaluate the action of the tidal torque in producing an angular momentum 
flux through wave excitation in section 5. In section 6 we describe numerical results while in section 7 we summarise the 
conclusions. 


2 BASIC EQUATIIONS 


We consider a uniformly rotating massive secondary star with mass Ms and radius Ra- We assume the angular velocity of 
rotation fig to be much smaller than the break up speed, i.e. ^ 1 , with = GMs/R^, so that effects of centrifugal 

distortion (oc 0,1) may be neglected in first approximation. 

We use spherical coordinates (r, 0, 93 ), with origin at the secondary’s centre, whereby 6 = 0 corresponds to its rotation axis 
which we assume to be parallel to the orbital angular momentum vector. However, we take the coordinates to be rotating 
with the secondary. 

The linearized hydrodynamic equations governing the non-adiabatic response of the uniformly rotating star to the perturbing 
potential ‘hr are the equation of motion 

^ + 2n»kxv' = —VP'+4 vP-V$t, (1) 

Ot p 

and the equation of continuity 

^+V-(pv')=0 (2) 

Here the velocity perturbation is v' and the pressure and density perturbations are P' and p' respectively. The unit vector 
along the rotation axis is denoted by k. 

We consider a close binary system in which the orbit is circular with angular velocity w and orbital separation D. The 
dominant tidal term of the primary’s perturbing potential is then given by the real part of: 

= fP Pi in) exp{2i[u}ft - p]) = (E>to exp{2i[u}ft - p]) (3) 


where Mp is the companion’s mass, cuf = uj — Os is the relative orbital frequency, as seen in the rotating frame, p, = cos9, 
Pi in) is the associated Legendre polynomial for I = \m\ = 2 and 

GMp 

■' 4P3 ■ 

For simplicity we have adopted the Cowling approximation, i.e. we have neglected perturbations to the secondary’s gravita¬ 
tional potential caused by tidal distortion. This approximation is reasonable because of the high central condensation. 

In this paper we adopt a simple model for calculating the tidal response in which the angular momentum exchange between 
the rotating star and the binary orbit occurs through angular momentum carying waves that are excited in the region of 
the convective core boundary and the radiative layers in its neighbourhood. It is supposed that these waves are subsequently 
dissipated with the consequent angular momentum exchange (see Papaloizou and Savonije 1984,1985, and Goldreich and 
Nicholson 1989). In order to analyse this model, only the interior regions of the star need to be considered. Then the 
perturbations are (almost) adiabatic such that 


dP' 

~dr 


TP dp' 
p dt 


+ v' 



(4) 


For simplicity in what follows below we shall take F to be constant and equal to 5/3. Then (^ and (^ can be combined to 
give 


dP' 

~dr 




(5) 


2.1 Reduction of the Response Equations 


The linear response to the tidal potential (^) is such that all perturbations have a p and t dependence through a factor 
exp(2i[cu/t — ip\). But note that, due to the occurrence of the Coriolis terms in the equations, the solution is no longer 
separable in 6, as it is for non-rotating stars. Equations (^) and (^) give for the components of the displacement ^ = (i(T)“^v/ 


- - 2i0sO sin = -K^-{^r + /g)A, 

or 

2 K dW 

- - 2iOsa cos9^,p = -—, 

r 09 


( 6 ) 

(7) 
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— + 2 iOscr (cos 61^0 + sin 6 '^r) 


K dW 
rsinO dip 


( 8 ) 


Here W = {P' + p^t)/P^^^, and K = P^^^ jp. The square of the Brunt-Vaisala frequency is given by 

VP /Vp_VP\ 
p V P TP J ’ 

the local acceleration due to gravity is given by p = —p~^{dP/dr), and the forcing frequency a = 2a;/. We can use equations 
(|) - (|) to express ^ in terms of W. Using the notation = {^ 1 ,^ 2 ,^ 3 ), (hi,/i 2 ,h 3 ) = (1,1/r, l/(r sin0)),and 

(* 1 , * 2 , * 3 ) = {r,8,p), we may write (using the summation convention) 




A 



dW , 
dxi g 



(9) 


The components of the Hermitian ( for real a ) matrix [H] are given by. 

All =+ 4a^Qs cos^ 9, H 12 = —sin 0 cos A 13 = 2 iQscr^ sind 

A22 = + A) + sin^ 9 , A23 = — 2 ifls(J cos 9 {—a^ + A), A33 = + A), 


with the unspecified components being given by the Hermitian condition assuming a to be real and 

A = + -^)(o'"* — 4cr^f2s cos^ 9) + 4a*Q,^ sin^ 9. 

Using the above expression for the components of the displacement in (P) gives a single second order partial differential 
equation for W in the form 


pi/r^ _ ^ _rp(i-i/r)^ 9 

qi dXi 



dW , 4 >t , 
Khj — -h Sji —A 

dxj g 


( 10 ) 


where (gi,®,®) = (r^,sin61,1). 

As the p and t dependence of W is through a separable factor exp{im{ujft — p)), with m = 2, after the replacement ^ —> —im 
and <I>T ^ 4?to, equation ( p^ becomes a second order partial differential equation of mixed type for W as a function of r 
and 9. It is hyperbolic whenever a real wave vector (fci, k 2 , 0) exists such that Aijkikj = 0. This condition is equivalent to the 
requirement that a satisfies the local dispersion relation (see Tassoul 1978, SPA) 

2 + 4Hs(fci cos^ — fe 2 sin6')^ 


(kl+kl) 


( 11 ) 


When this can be satisfied and the boundary conditions are non-dissipative, one expects a dense spectrum of normal modes 
and a singular response to forcing (SPA). This problem is avoided here by, from now on, allowing a to have a small negative 
imaginary part. This latter Landau prescription corresponds to the forcing potential being slowly switched on at time t = — 00 . 
This also leads naturally to the selection of a predominantly outgoing wave boundary condition, corresponding to the physical 
situation where the surface regions of the star are assumed to be highly dissipative with little or no wave reflection occurring 
from them. 


3 ASYMPTOTIC TREATMENT FOR LOW FREQUENCIES 

Because of lack of separability, the solution of equation (^^ cannot be undertaken analytically in general. However, progress 
can be made if we assume that the square of the forcing frequency is small compared to the Gp, the square of the inverse 
of the local dynamical timescale. We shall be specially interested in the inertial regime and so we shall assume is 

comparable to unity. 


3.1 The radiative zone 


In a radiative region A is non zero and positive and for low frequencies jA can be considered to be small, apart possibly 
from a region in the vicinity of the interface between the convective core and surrounding radiative zone. However, what 
follows below can be shown to be valid also in such a region. We write equation (^^ in the form 


pi/rp^_^rp(i 


qi dxi 




= s. 


( 12 ) 


where we have grouped the forcing terms involving <E>to together on the right hand side of (|l2| ) in S. Considering that e = j A 
defines a small parameter, we may attempt to find solutions by expanding the left hand side of (^|) to first order in this 
parameter. We also allow the solution to vary rapidly with radius such that ^ = 0(e“^^^)i. It is further assumed that the 
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angular variation is much less than the radial variation which is the situation for the expected excited low frequency g modes. 
Then expansion of ( p^ ) in powers of e gives to zero order 

^ 2 pp(i-i/r) Q f dW\ KVP , 


P^^^WP - 

dr y A dr 

where the operator 0± is defined through 

„ /Ti/N _ ‘^mQsfJ.Q m'^W 

^ dg a{l — g'^) — 


+ 


-Ox (IT) = PS, 


with Q being related to W through 

2,dW 2mQsgW 


Q 1- 


40 ? m " 


= -(1 - M 


dg 


(13) 


(14) 


(15) 


3.2 The convective core 


We idealize the convective core to be a region in which .4 = 0 interior to the boundary, where we suppose .4=0, but V.4 7 ^ 0. 
When .4 = 0, we cannot make an expansion based on the smallness of cP jA. In this case we must retain (|^). However, S 
takes the simple form S = p 4 >T 0 , and thus (p^ becomes 


pi/r^ P — 

qi dxi 




(16) 


Noting that .4 = 0, the unforced form of ( p^ ) gives rotationally governed inertial modes in the core (see Papaloizou and 
Pringle 1981, SPA). The local dispersion relation is ( pi] ) with .4 = 0. In order to solve (@) boundary conditions for W are 
needed. In the low frequency limit we may use the fact that the stratification in the radiative zone leads to being small 
compared to the other components of the displacement. This is a familiar property of oscillation modes such as g modes in 
radiative layers. This suggests the boundary condition = 0 at the convective core boundary. Using (^, this can be expressed 
as a condition on W. 

For real a finding W in the convective core is problematic because of the dense spectrum of the inertial modes (Greenspan 
1968, Papaloizou and Pringle 1981, SPA). However, if a is complex as we shall assume, singular responses are avoided except 
possibly when there are near resonances with global modes. A similar situation is expected if viscosity is introduced. Solving 
(@) subject to = 0 on the convective core boundary is equivalent to solving an elliptic boundary value problem and the 
solution will provide W on the core boundary which may be used as a boundary condition there for the problem of finding 
IV in the radiative zone. 

The elliptic boundary value problem described above is not simply soluble in general, even in the low frequency limit. However, 
an exception occurs when the core is assumed to have constant density and pressure, an assumption which is only justifiable 
when the core is small. Then a solution can be found in the form of a finite polynomial. To find such solutions, it is best to 
write equation in cylindrical polar coordinates (w, ig, z) in the form 


pa 


TP 

where 


W + - 


1 

■ 1 d 


rrdW' 

d^w 

{l-x^) 

■uj dvj 

V dw) 

■w‘^ 

dz^ 


rp(i+i/r) ‘ 


(17) 




In order to proceed we take P and p to be constant and equal to their values on the convective core boundary. In the low 
frequency limit the first term on the left hand side of (|l^ may be neglected. The problem is then equivalent to calculating 
the response of an homogeneous incompressible sphere (see Greenspan, 1968). In this case the solution that has zero normal 
displacement at the core boundary can be written as a polynomial in vj and 2 for m = 2 in the form 


Ti/ 4,7 2 I 22 

W = actZJ -|- OcU7 + CcZ VJ , 


where 


fcc = - 


(18) 


3fp^a^ {l — x^){2 + x) 
= 2 rp(i+i/r) (14 -x^ + 7x) ’ 

SfpGla^ {l + x)G-x^) 

2 rP(i+i/r) (i4_a;2 + 72;) ’ 

3fp^a^ ( 1 -I-a ;)(2 — 2;) 

= 2rP(i+i/r) (14 _ 2;2 72;) ’ 


and 
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where is the radius of the convective core. We remark that the response is resonant when 14 — + 72; = 0. There are 

thus only two real resonant forcing frequencies corresponding to the roots of the quadratic given by x = (7 db (105))/2, or 
Ljj/Q,s = (±-^(105) — 7)/28. Although the spectrum is dense, the nature of the forcing potential allows only two potential 
resonances to be excited in this case. It would seem reasonable to assume that other modes can only be weakly excited in 
more general cases. In fact the above solution for a uniform core may be regarded as the first member of a sequence of higher 
order polynomial approximations to solutions in more general cases (see Papaloizou and Pringle, 1981). Accordingly we shall 
use to determine W on the core surface in the analysis presented below. 


4 ASYMPTOTIC SOLUTION IN THE RADIATIVE ZONE 

The governing equation (0) developed for small jA can be solved by utilising the fact that the operator on the left hand 
side is separable (see SPA). We may write 


IU = ^w„(r)X„(Ai). (19) 

n=0 


Here the A„ are a set of orthogonal functions of /r over (—1,1) so that the expansion coefficients are given by 


WXnd^ 



The Xn satisfy a second order ordinary differential equation defined through the first order pair (see SPA) 


dV„ , mxfiVn m^Xn _ ^ ^ 

~ (1-^2) ” (1 -M^) “ 

Vn (l - = -(1 - - mXflXn- 


( 20 ) 

( 21 ) 


The eigenvalues An are determined so as to make regular at — XI. Note that they are functions of the parameter x and 
thus the forcing frequency and the stellar rotation rate. Note too that for potentials of the type we consider here which are 
even functions of fi, we may restrict ourselves to An such that Xn,T>n are even and odd functions of /r respectively. 

Using (^ we obtain 


0±{W) = -J2 AnU;„(r)Xn(M). (22) 

n=0 

Using the expansion given by equation (^^ and the orthogonality of the X^, equation (§) gives 

(T^pr^\ d f r'^KP^/^ dwn\ _ a'^pr^Sn 

VP J KPXr ar A dr J ~ ppi+i/r ’ ^ 

where 

fh SX^dp 

Sn = -Tj . 

J_,X^dp 

The problem is thus reduced to solving a set of second order ordinary differential equations for Wn{r). This is possible because 
of the fact that the normal mode problem becomes separable in the traditional approximation (see Chapman and Lindzen, 
1970 for a discusion in the context of atmospheric tides and also SPA). This approximation neglects the 6 component of the 
stellar angular velocity (Unno et al 1989), and it is expected to become valid in the stratified radiative layers in the limit of 
low ferquencies when the radial displacement becomes small compared to the angular displacements. 



4.1 Reduction of S 

From direct calculation using equation (0) we obtain for large A 

ppi-i/r d /_ x^VP d f p{l - p‘^)^To \ rraxPR-Fro 


dr 


gr dp Y — J gr(l — x'^p^) 


+ P^TO , 


and after performing an integration by parts 


Xidp= / SXr,dp = fpr" 


d x^VP 




dr 


■ 


pgr 


J + T 


(24) 


(25) 
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with 


and 


1= J Xr,Pi{n)dfl, 



mXn/x)Pi{iJ,)d^. 


4.2 Behaviour of the eigenvalues Xn{x) 

The calculation of the tidal torque has thus been reduced to solving the sets of second order differential equations The 
boundary conditions are that the solution should, if possible, correspond to a predominantly outgoing wave at large radii, 
and Wn should be matched at the core boundary to the solution obtained above that applies to the convective core. 

Clearly the form of the solutions of the determining equation depends on the eigenvalues An- /.From general scaling 
arguments one expects Wn = 0{P/Xn) which vanishes as a vanishes. When An ^ 0, the unforced solutions for Wn are short 
wavelength g mode like waves. There is special interest in the smallest positive An because these lead to solutions with the 
longest radial wavelengths which will give the strongest responses to global forcing of the type considered here. 

For a non-rotating star {x = 0) the eigenvalues X„ = l{l + 1), with I = n = 2,4,6... and the corresponding eigenfunctions 
Xn — P^{fj.). However, the behaviour of X„ is very different in the inertial regime for which —1 < x~^ < 1. We plot the 
smallest positive An we obtained for a range of values of x in the inertial regime in figure 1. It will be seen that there is a 
significant difference between positive and negative x. For a; > 0, corresponding to prograde relative rotation of the companion. 
An is similar to the non-rotating case. However, for a; < 0, corresponding to retrograde relative rotation, An becomes large 
for x~^ < —1/6. These eigenvalues correspond to disturbances conhned cosely to the equator. For —1/6 < x~^ < 0, another 
small An exists corresponding to a toroidal mode resonance. These eigenvalues are important because they lead to solutions 
with strong global responses to the forcing tide. 


4.3 Toroidal mode resonances 

Strict toroidal mode resonance can be defined to occur when x is such that an eigenvalue An = 0, (see SPA). Then we have 
X = —{l{l + l))/m, with I = n > \m\ being an odd integer to provide solutions with the required symmetry type. These 
resonances corespond to the r mode frequencies (Papaloizou and Pringle, 1978) 

—2mfls 

The corresponding eigenfunctions are T>n = P;'"*'(/i), and X„ = m~^(—l{l + l)g,Pl^\g,) — {1 — g,^)dPl”^^ {g.)/dg,), where P(™''(/r) 
denotes the standard Legendre function. The lowest order resonance has I — 3 so that x~^ — —1/6. 

We see that, for positive values of m , the toroidal mode resonances occur for negative forcing frequencies, corresponding to 
retrograde orbital rotation relative to the star. In this case the star rotates faster than the orbit so that the action of the 
tides is to cause the star to spin down. It is clear from the above discussion that the existence of toroidal r mode resonances 
leads to a qualitatively different tidal response for some negative forcing frequencies which is in general much larger than for 
corresponding positive forcing frequencies. The tidal evolution timescale will accordingly be shorter in that case. 


5 WKB SOLUTION 

It is possible to solve equation (^^ by a Green’s function method as in Papaloizou and Savonije (1985). In general extended 
regions outside the convective core are found to contribute to the excitation of an outgoing wave. However, for simplicity we 
here look at the solution of equation (^^ in the low frequency limit where a WKB approximation should be adequate. In this 
case the largest lenghthscale associated with the solution for occurs near the convective core boundary (Zahn, 1977) where, 
for An > 0, an outgoing wave is excited. Only this region matters for wave excitation in the limit of low forcing frequencies 
but if this limit is to be of practical use, the convective core should not be too small. 

To construct the solution, we write z = r — r^, with being the radius of the convective core boundary. Further we assume 
a local first order Taylor expansion such that A = {A!)cZ, and evaluate all other quantities in (|^ at the core boundary. We 
then obtain 
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Here the subscript c denotes evaluation at the core boundary. From now on we shall take this as read for all state variables 
and the subscript c will be dropped. 

It is straightforward to express the solution of ( p^ ) that coresponds to only outward going waves in terms of Hankel functions, 
H^\ as follows. If Q = z~^dw„/dz, then the solution for Q is given by 


35/6r(4/3) exp(-j 7 r/ 6 ) ., 1/2 ^( 1 ) /2 i/ 2 ., 3 / 2 \ 

Q 2 ai /6 ^ 1/3 [ 3 °^ ^ J’ 

(27) 

where 


A' ^ a^pr^] 

0 - 2^2 " TP J ’ 

(28) 

Co =-arun(O) - 

(29) 


with ic„( 0 ) denoting the prescribed value of on the core boundary. 

Here it is assumed that the real part of the forcing frequency is positive. Solutions appropriate to negative forcing frequencies 
can be found by setting m = 2 —> —m, below while retaining the forcing frequency as positive. 

The Wn{0) are given by 


w„{ 0 ) = 


WXud^ 


with W being evaluated on the core boundary. Using the solution given by (p^, we obtain 

+ bcr^(l - fj.^) + Ccr*fi^{l - iJ?))Xnd^. 


w„{0) = 




5.1 Angular momentum flux and tidal torque 


The asymptotic outgoing wave solutions outlined above are associated with a conserved angular momentum flow or wave 
action ( Goldreich an Nicholson, 1989). Assuming the outgoing waves are ultimately absorbed near the surface, this angular 
momentum is deposited there and correspondingly removed from the orbit. But note that this exchange can be of negative 
sign for retrograde forcing. Thus in this picture, the tidal torque acts through the production and absorbtion of angular 
momentum carrying waves. 

The radial component of the wave angular momentum flux appropriate for responses to the complex forcing potential is given 
by 

F = jIM (P'C), 

where XA4 denotes the imaginary part (see Ryu and Goodman, 1992, Lin et al, 1993). Using 

Cr = -and P' = 

in the asymptotic limit, we find the asymptotic form ( for large 2 ; ) of the angular momentum flux associated with each An 


F = 


38/3 (r(4/3))^ 

87rai/3(^')2 (^An- 


(30) 


The total rate of change of angular momentum of the star is found by integrating the radial flux over a spherical surface at 
the convective core boundary. It is thus given by 

/•i 


J = 27rr'^ 


Fdfi. 


Using this, while neglecting {a‘^pr^)l(TP) in comparison to An, we find after some algebraic manipulation that for a particular 
An, 


j = pr^tT^l 


3«/3(r(4/3))=e2|H|2 f 


27/3An/"/n 


He 


TP 


rA' 


( 31 ) 


Here 


MlRl 

MlD^ ’ 
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and 0 is defined through 


S„ = - 


fpP 


TP (AnpP 


pgr 


+ J 


pgr 


fpP 


0 , 


with M being the mass interior to radius r and 


/n = 


Xldp. 


The quantity ^ is defined by 


H = 1 + 


w„(o)A„rpi+irj^ 

fp2rifj2Q 


Here we have set m = 2, corresponding to prograde forcing frequencies. In this case J > 0, corresponding to stellar spin up. 
For negative forcing frequencies the same expression may be used but the sign of J is reversed. In general we should sum 
over \n. However, for simplicity we shall restrict ourselves to the values displayed in figure 1 for the numerical work described 
below. 

In order to compare with the work of SP, we introduce to defined through 

kMsRi^T\u!s\ 


to 


with kMsRs being the stellar moment of inertia. Thus 


to 


kMsRi 




(\pA 

Fl^pr^ 38/3(r(4/3))^02|H|2 ^ 


-5/3 /n2 2 

^cpr 

VP 


J 


rA' 


- 1/3 


(32) 


We remark that for a particular stellar model 0^|H|^ = F{x) depends on the forcing and stellar rotation frequencies through 
the quantity x only. Thus we have the scaling law that 


tooc(|w/|) /F{x), 


(33) 


which, as long as our simple prescription for including damping at inertial core mode resonances is used ( see below ), enables 
results to be scaled to different forcing frequencies for the same value of x. 


6 NUMERICAL EVALUATION 

We have evaluated to using ( |3^ for the 20 Mq stellar model considered in SP. We display the result for logj^g to plotted 
against uif/Fla for Sis = O.lSlc in the inertial regime in figure 2. But note that these results may be scaled to other values of 
FIs using the scaling relation (|3^. 

Although we do not reproduce all the details found in SP, the most important features where they can be compared, including 
the form of the difference between prograde and retrograde forcing, and agreement with the values of to, to order of magnitude, 
are found. This type of qualitative agreement is reasonable in view of the large amount of variation in to and the fact that 
the asymptotic analysis is carried out in the limit of both low forcing and rotation frequencies. This results in extreme 
sensitivity to the location of the convective core boundary. Also simplifying approximations were made to find the response 
of the convective core. Note too that the assumption of an outward going wave condition in the low frequency limit has 
prevented the appearance of the oscillatory behaviour found in SP. It is hoped that a future extension of the analysis will 
enable relaxation of these approximations. 

We do find a strong toroidal mode resonance at the expected location u>f = —IIs/6. This occurs because of the long 
wavelength response near this location resulting from the small Xn. We remark that the assumption of an outgoing wave 
condition will break down in the centre of this resonance where the response is of very long wavelength. Also we find rather 
small values of to near x~^ = 0 ( note the resolution here was 0.01) in a domain that could not be considered in SP. The 
origin of this behaviour can be traced to the forcing term oc x'^ in (|25|). 

As in SP, we obtain larger to in general for w/ < —fIs/6 because of the larger A„. Note too that we have two core inertial mode 
resonance features at uif ^ —O.Ofls and w/ ~ 0.12IIs. Similar features were seen in SP. To obtain features of similar scale, 
we incorporated damping in the core by adding a negative imaginary part to w/ of magnitude 0.06|Re(t<;/)| when calculating 
uin(O). In reality the damping is likely to be due to g mode losses, that would appear in an improved higher order theory, 
as well as applied viscosity. With our simple damping prescription, the details of these particular features are not too well 
reproduced but this is not too surprising in view of the simplifying assumption of an homogeneous core. 
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7 CONCLUSIONS 

We have developed a generalization of the asymptotic treatment of the tidal torques experienced by a massive star as a result 
of a companion in circular orbit considered by Zahn (1977) and Papaloizou and Savonije (1985) to the case of a rotating star 
when the forcing frequency is less than the rotation frequency in magnitude. The results confirm the presence of a strong 
toroidal mode resonance feature for retrograde forcing and also, albeit with a very simplified model of the convective core, 
the presence of some core inertial mode features found in SP. 
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Figure 1. The form of logj^Q as a function of x~^ used in the evaluation of to - These An represent the lowest positive values outside 
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Figure 2. A plot of logj^Q to? ^0 y^-? as a function of x ^ in the inertial regime for fls = O.lflc. 
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